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We study the Josephson plasma resonance (JPR) in Bi2Sr2CaCu20s+i/ (BSCCO) with spatially 
dependent interlayer-phase coherence (IPC). The half-irradiated BSCCO (HI-BSCCO), in which 
columnar defects are introduced only in a half of the sample, shows several resonance peaks, which 
are not simple superposition of the peaks in irradiated- and pristine-parts. JPR in HI-BSCCO 
changes its character from irradiated- to pristine- type at a crossover frequency (ui C r)- We demon- 
strate that the one-dimensional linearized sine-Gordon equation, which takes into account the spatial 
dependence of IPC, can reproduce most of the experimental findings including the presence of ui cr - 

PACS numbers: PACS numbers: 74.60.Ec, 74.72.Hs, 74.25.Dw, 74.25.Ha 



Highly anisotropic high temperature superconductors 
like Bi 2 Sr 2 CaCu208 + j ; (BSCCO) can be considered as 
natural stacks of Josephson junctions. In such systems, 
collective oscillation of superfiuid density gives rise to 
the plasma oscillation along the stacking direction of o 
axis . A resonance occurs when the frequency of the 
incident electromagnetic field matches the plasma fre- 
quency. This resonance is called Josephson plasma res- 
onance (JPR) and its frequency is give by the following 
equation 

Up = w p0 < cos 4>n,n+l > ■ (1) 

Here 4> n ,n+\ is the gauge-invariant phase difference 
between the neighboring layers which determines the 
interlayer-phase coherence (IPC), lu p q (= c/\ c ^/e$) is the 
Josephson plasma frequency at zero field (c, A c , and eo 
are light velocity, c-axis penetration depth, and dielec- 
tric constant, respectively) and < > denotes thermal 
and disorder average . JPR gives valuable information 
on the charge dynamics, superconducting symmetry, and 
vortex states |j, la, Ifl 151 la ISl lill m anisotropic super- 
conductors. 

Most of JPR experiments, however, have been per- 
formed in macroscopically homogeneous systems. What 
would happen when we have different coupling strengths 
in a system? In this context, JPR in T*-phase super- 
conductors, which have two Josephson junctions with 
different coupling strengths in the unit cell, has been 
studied extensively both experimentally and theoretically 
[Til H2I. Ik3. IbH llj| . Experimentally found two plasma 
modes were analyzed in the framework of Josephson junc- 
tion stacks with alternating coupling strengths. Another 
interesting question is what would happen if IPC has 
in-plane distribution on a macroscopic scale. JPR fre- 
quency is greatly modified in the presence of vortices, 
since the configuration of pancake vortices determines 
< cos0 n! „+i > in Eq. (1). When the pancake vortices 



are aligned along the c-axis, < cos0„ jrl +i > is equal to 
unity and lo p takes the maximum value. A strong sup- 
pression of < cos0 n! „+i >, and hence of u> p , occurs in 
magnetic fields due to the misalignment of pancake vor- 
tices. Columnar defects (CDs), generated by heavy- ion 
irradiation, act as directional pinning centers and dra- 
matically enhance < cos</>„ iI1+ i > 0,11 H- We can cre- 
ate an arbitrary distribution of IPC by controlling the 
pattern of heavy-ion irradiation. In this paper, we re- 
port JPR in the simplest inhomogeneous system having 
two regions with and without CDs. We find complicated 
field dependence of resonance peaks including their num- 
bers. We present a clear explanation of this non-trivial 
behavior of JPR through numerical calculations of one- 
dimensional sine-Gordon equation. 

Single crystals of BSCCO are grown by using 
the traveling-solvent-floating-zone method hj. The 
optimally-doped crystals (T c = 90.0 K) are irradiated by 
600 MeV iodine ions (A) (JAERI) and 6 GeV Pb ions 
(B, C and D) (GANIL) at the matching field of 5$ = 
20 kG. The thickness of these crystals are small enough 
so that the heavy-ions can pass though. The BSCCO 
crystal is half-irradiated (HI-BSCCO) by covering a half 
of the crystal with a gold foil, which is thick enough to 
stop the incident heavy-ions into the sample. Dimensions 
of the samples are 670x670x20 ^m 3 (A), 670x670x20 
^m 3 (B), and 1480x670x20 ^m 3 (C, D). Samples C and 
D have the same dimensions but different distributions of 
CDs (see Fig. 3). To study sample size dependence, sam- 
ple B, C, and D are irradiated at the same dose. JPR 
is measured by the cavity perturbation method at mi- 
crowave frequencies of 24.2, 41.7, 46.3, 49.4, 52.0, 54.9, 
56.6, and 61.2 GHz |y. The crystal is set in a cylindrical 
copper cavity so that the microwave electric field is ap- 
plied parallel to the c-axis. The external magnetic field 
up to 90 kOe is also applied parallel to the c-axis by a 
superconducting magnet. 
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Figure 1 (a) shows microwave power absorption as a 
function of field, P abs {H), for HI-BSCCO (A) at 50 K for 
various frequencies. The dashed lines in Fig. 1(a) show 
JPR fields in irradiated and unirradiated samples with 
the same doping levels. JPR in HFBSCCO shows a com- 
plicated behavior, which is not a simple superposition 
of resonances from the constituent parts of the sample. 
Change of frequency does not only shift the resonance 
field Hp, but also changes the number of resonance peaks. 
At 24.2 GHz, only a single resonance peak appears close 
to the H p in the irradiated sample. As frequency is in- 
creased from 24.2 GHz to 61.2 GHz, H p of the main peak 
shows characteristic anticyclotronic behavior [5J. Small 
shoulders observed at 41.7 and 46.3 GHz at the low field 
side of the main peak are most probably caused by orig- 
inal inhomogeneity of the crystal |17| . Above a crossover 
frequency of w cr /27r ~ 55 GHz, H p suddenly shifts to 
lower field. It would be worth mentioning that when we 
measure JPR for a pristine and an irradiated (_B$ = 20 
kG) samples placed side by side, the observed resonance 
peaks are simple superposition of the peaks in each indi- 
vidual sample (not shown). This indicates that the elec- 
tromagnetic continuity condition between pristine- and 
irradiated-part plays an important role to the behavior 
of JPR in HI-BSCCO. 

Spatial and temporal dependence of the phase in a 
layered superconductor is described by a coupled sine- 
Gordon equation^]. According to Koshelev |19j. when 
the IPC along c-axis is uniform, spatial variation of the 
amplitude of plasma oscillation (8(x)) obeys the lin- 
earized sine-Gordon equation (LSGE) 

(^^-c( x) )e( x) + xtvle( x ) = -^ (2) 

where u) and C{x) represent incident microwave fre- 
quency and spatial dependence of IPC, respectively. v Cl 
D, and Jo (= < I > oWpo e o/87r 2 s) are damping frequency, mi- 
crowave amplitude, and maximum Josephson current at 
zero field (<I>o and s are the flux quantum and the inter- 
layer spacing). In the present case, the spatial distribu- 
tion of C(x) is 

= f Ci < x < L x 

\C 2 L 1 <x<L 1 +L 2 (=L)' 

Here, L\ and L 2 are widths of pristine and irradi- 
ated parts (see inset of Fig. 1(b)) and L is the total 
width of the sample. Boundary conditions for Eq. J5J, 
6>'(0) = 6'(L) — 0, are the requirement from no bias 
in-plane current at sample edges. Microwave absorption 
power is originated from Joule-heating, and is given by 
P{oj) oc Sj- Jq |6'(a;)| 2 da;. In the case of a homogeneous 
system where C (x) in Eq. J3J) has a fixed value C all over 
the sample, the resonance occurs at a single frequency of 
u)p = CcjpQ. Figure 1(b) shows the calculated absorption 
power as a function of external field. In this calculation, 
field dependence of C(x) for each part is assumed based 



on our previous results |6J. It should be noted, however, 
since the exact field dependence of IPC in the irradiated 
region is not known, we approximate it by three straight 
lines in the double logarithmic plot as shown by dotted 
lines in Fig. 2 [2fJ . We also assume appropriate values of 
v c , tOpo, cq and A c as shown in Fig. 1(b). A single reso- 
nance peak appears near the resonance field for the irra- 
diated region below 45 GHz and the resonance peak shifts 
rapidly towards lower field between 50 GHz to 60 GHz, 
which we identify as w cr /27r. The LSGE reproduces all 
characteristic features of JPR in HI-BSCCO well. Slight 
quantitative disagreements could be caused by the non- 
monotonic field dependence of IPC in the irradiated part 
at low fields |3,|2J. On the other hand, above uj cr , new 
weak peaks start to emerge at higher fields. These peaks 
appear since Eq. (2) without the damping term has the 
same form as the Schrodinger equation. The intensity 
of the new peaks becomes stronger at higher frequencies. 
Since our measuring frequency is limited, we cannot ob- 
serve the new peak in this sample. 

Figure 2 shows calculated field dependence of uj v in 
HI-BSCCO with the sample size of L = 670 /im (sample 
A) and L = 1480 /im (sample D). We also add experi- 
mental data as solid circles in Fig. 2. In this calculation, 
HI-BSCCO with L = 670 /im shows three peaks in the 
frequency range, < (cu/ujpo) 2 < 1. In the case of iso- 
lated system, resonance peak in irradiated and the pris- 
tine parts are observed on C\ and C 2 line, respectively. 
Within our microwave frequency window, only the lowest 
frequency resonance is observed for L — 670 /im sample. 
Still, calculated results agree quantitatively with experi- 
mental data. The calculated field dependence of the uj p 
for HI-BSCCO with L = 1480 /im is also shown with red 
lines. Since doping levels in sample A and D arc different, 
C 2 in sample A and D are shown as blue and red dotted 
lines, respectively. To avoid complexity, we draw only 
two red lines as the lowest and the next lowest uj p for L 
= 1480 /mi sample. In addition, it is clear that oj cr /27T 
is decreased to ~ 30 GHz with increasing the sample size 
perpendicular to the boundary. 

To confirm the size dependence of oj cr , we measure 
sample (B)-(D) with the same doping levels to each other. 
Here C(and D), which have twice the size parallel (and 
perpendicular) to the irradiation boundary of sample (B), 
respectively. Figure 3 shows P a bs{H) in these samples for 
41.7 GHz at 50 K. The resonance in sample C occurs at 
almost the same field as in sample B. This fact is valu- 
able in further convincing ourselves that it is sufficient to 
consider only the one-dimensional equation as Eq. (2). 
The resonance peak in sample D, which corresponds to 
the above calculation with L = 1480 /im, shifts to lower 
field at 2.5 kOe and a new peak appear at 26 kOe. In 
Fig. 2, we plot these peaks as closed squares. The lower 
field peak and the new peak locate favorably on the low- 
est Up and the second lowest oj p curve. These results 
strongly suggest oj cr /2iT in sample D is reduced below 
41.7 GHz region and support our claim that uj cr tends to 
decrease with increasing the sample size perpendicular to 
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the boundary. 

Let us now discuss the behavior of JPR in HI-BSCCO 
in more detail. For this purpose, we calculate the lu p as 
a function of C'i with fixed C\. Figure 4(a) shows the u> p 
as a function of C 2 X (0.001 < C 2 < 1) with fixed d = 
0.001. Since IPC changes more significantly in the irra- 
diated part in our experiments, this figure approximates 
the actual field dependence of u p . Dotted lines in Fig. 
4(a) represent the conditions of uj 2 /oj 2 = C\ and Ci. At 
the intersection of the two lines (1), namely C\ = C 2 = 
0.001, HI-BSCCO becomes uniform and shows only one 
resonance at oj 2 = Ciiv p0 = 0.001Wp . The increase of Ci 
causes an increase of iv p (•), and in addition, brings two 
new resonance peaks at higher frequencies (■) and (A). 

In order to get further insight into the nature of the 
resonances, we calculate spatial variation of the real part 
of 9{x) (Re[0(x)]). The Eq. (2) is now analogous to the 
Schrodinger equation for a particle in spatially dependent 
potential well with modified boundary conditions. The 
Re[0(x)] shows a snapshot of the resonance state in the 
sample. Figure 4(b) shows Re[6*(x)] normalized by its 
maximum value in HI-BSCCO at several values of C 2 l 
and oj p ((1) to (6) in Fig. 4(a)). Re[0(x)] is constant for 
homogeneous case of (1) as expected, whereas it shows 
spatial variation for C\ 7^ Ci as in (2)-(6). When C 2 is 
slightly larger than C\ as in (3) and (6), Re[0(x)] changes 
sinusoidally and is antisymmetric (or symmetric) with 
respect to the center of the sample. As Ci increases fur- 
ther, sinusoidal variation of Re[#(x)] is modified. When 
iVp/iVpQ is equal to C2 as in (4), Re[0(x)] changes sinu- 
soidally only in the pristine part and is nearly constant 
in the irradiated part, since the wave length of Re[0(x)] 
in the irradiated part diverges to infinity. In the case of 
C\ < uj p /uj p0 < C2, the spatial variation of Re[0(x)] in 
the irradiated part is damped as shown in (2) and (5). 
The finite amplitude of Re[0(x)] in the irradiated part is 
the result of the continuity condition at the boundary. 

To clarify the origin of additional resonance, we calcu- 
late the lo p ignoring the damping term in Eq. (2), which 
has little influence on the determination of lo p . The re- 
sulting equation has the same form as the Schrodinger 
equation for a particle in a box. In an uniform system 
(Ci = C2), spatial derivative of is neglected, and we 
have only one resonance at Lu 2 /tu 2 = C\ = Ci. In case 
of u> p /u) p0 > C\ 7^ C2, 0(x) is represented by the sum of 
transverse- and longitudinal-components as 

with (j=l, 2). Here kj (= (1/A C ) ^Juj 2 /^ - Cj) is the 

wave number. In the case of C\ ~ Ci, 9{x) is approxi- 
mately expressed by a sinusoidal wave as in (3) and (6), 
since k\ and ki are nealy equal. Then k\ and ki must 
satisfy the condition sin(fciL) <~ sin(fc 2 L) = by the 
boundary condition of Eq. (2), from which we can eval- 
uate uj p as u> n , where 

^njCliL) A c 2 

-2 = Ci + (— mr) (3) 



with n = 0, 1, 2, • • •. In the case of Ci > ui 2 /u> 2 > C\, 
lu p can be obtained from the intersection of two lines 
ki tan(fci_L/2) = k 2 tanh(« 2 i/2) and n\ + k\ = C 2 - C\ . 

Here n 2 = (l/X c )^Jc 2 — uj 2 /w^q is the attenuation coef- 
ficient in the irradiated part. In the limit of C 2 3> C\, 
the two lines cross to each other at k\L/2 — ir/2 + mr, 
which translates into k\ = ir(2n + 1)/L. This means 
that uj p started from ui n (n = 0, 1, 2, • • •) approaches to 
u)i n +i(Ci,L). Dotted lines in Fig. 4(a) indicate ui p — 
w n {Ci,L) for n — 1,2,3 with L — 670 /an and C\ = 
0.001. Actually, in the limit of Ci — > C\, each resonance 
frequencies (•, ■. and ▲) approaches to w n with n = 0, 1, 
and 2. By increasing Ci, u> p started from Lo n shift towards 
higher frequency and it approaches to a> 2 n+i in the limit 
of large Ci. It should be noted here that in the real 
system, the maximum of Ci is limited by unity. There- 
fore, JPR frequencies in HI-BSCCO can be classified by 
characteristic frequencies to n {C\,L)(n = 0, 1, 2, • • •). 

Finally, we discuss the origin of uj cr . The observed 
plasma resonance in HI-BSCCO with L — 670 /im cor- 
responds to the lowest uj p (•) in Fig. 4(a), since we mea- 
sure at relatively low frequency In the calculation with 
fixed C\ shown in Fig. 4(a), the lowest oj p approaches to 
u>i(Ci,L) in the limit of Ci ^ C\ at lower fields. How- 
ever in reality, the limiting value u>i itself is a function 
Ci(cx l/H). Considering the field dependence of C\, we 
calculate uj 2 (Ci(H), L)/lu 2 for L = 670 and 1480 jitm, 
which are shown in Fig. 2 by dash-dotted lines. Now it is 
clear that the crossover frequency uj cr appears since uj 2 
follows ijj 2 (Ci(H), L) at lower fields when Ci is almost 
saturated. On increasing field, uu p (H) starts to follow 
the lowest line in Fig. 4 since the change in C\ is much 
less than that in Ci. 

In summary, we study the Josephson plasma reso- 
nance in half-irradiated Bi2Sr2CaCu208+ 2/ (HI-BSCCO) 
with inhomogeneous phase coherence. All experimental 
results, including multiple resonance and the crossover 
frequency, are well reproduced by calculations based on 
the linearized sine-Gordon equation. Resonance peaks in 
HI-BSCCO are classified in terms of characteristic fre- 
quencies, which belong to different modes with specific 
wavelength relative to the sample size. In such layered su- 
perconductors with spatially dependent interlayer-phasc 
coherence, microwave electric field parallel to the c-axis 
can partly excite the transverse mode, which can only 
be excited by microwave magnetic field in homogeneous 
systems. 
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FIG. 1: (a) P abs (H) for HI-BSCCO at 24.2, 41.7, 46.3, 49.4, 
52.0, 54.9, 56.6, and 61.2 GHz at 50 K. Dashed and dotted 
lines indicate the resonance fields for pristine and irradiated 
part of the sample, respectively, (b) P a bs(H) calculated based 
on the linearized sine-Gordon equation (Eq. (2)) for position 
dependent IPC. 
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Fig. 2 



FIG. 2: (Color online) Field dependence of uj p for HI-BSCCO 
with L= 670 (o) and 1480 fim (red lines) calculated by lin- 
earized sine-Gordon equation. Solid circles (•) and Solid 
squares (■) show observed resonance peaks at 50 K. The 
dashed line and dotted lines show C\(H) and C2(H) which 
are used in this calculation, respectively. For explanation of 
LO cr and lu\{C\,L), see text and Eq. (3). 
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Fig. 3 



FIG. 3: P aba (H) in three HI-BSCCO samples with different 
dimensions and irradiation patterns with B$ = 20 kG for 41.7 
GHz at 50 K. Resonance peaks are marked by arrows. The 
irradiation pattern for each sample is schematically shown. 
The difference between samples C and D is the direction of 
the irradiation boundary. 




Fig. 4 

FIG. 4: (a) C2 dependence of the peak frequency (ui p ) (solid 
symbols). Ci is fixed at 0.001. Dotted lines represent Ci and 
C2. Dashed lines represent uj^/ui^q (n = 0, 1, 2, and 3) with 
Ci = 0.001 and L\ + L2 = 670 /im. Here, ui n — u n (G\,L). 
(b) The distribution of Re[0(a;)] in sample at various points 
in (a). 



